SARS-CoV-2 spike variants differ in their allosteric response to linoleic acid.

The SARS-CoV-2 spike protein contains a functionally important fatty acid (FA) binding site, also found in some other coronaviruses (e.g. SARS-CoV and MERS-CoV), which binds linoleic acid. When occupied by linoleic acid, it reduces infectivity, by 'locking' the spike in a less infectious conformation. Here, we use dynamical-nonequilibrium molecular dynamics (D-NEMD) simulations to compare the response of spike variants to linoleic acid removal. D-NEMD simulations show that the FA site is coupled to other, some distant, functional regions of the protein, e.g. the receptor-binding motif, N-terminal domain, furin cleavage site, and regions surrounding the fusion peptide. D-NEMD simulations also identify the allosteric networks connecting the FA site to the functional regions. Comparison of the response of the wild-type spike with four variants (Alpha, Delta, Delta plus, and Omicron BA.1) shows that the variants differ significantly in their response to linoleic acid removal. The allosteric connections to the FA site on Alpha are generally similar to those on the wild-type protein, with the exception of the receptor-binding motif and the S71-R78 region, which show a weaker link to the FA site. In contrast, Omicron is the most affected variant exhibiting significant differences in the receptor-binding motif, N-terminal domain, V622-L629, and the furin cleavage site. These differences in allosteric modulation may be of functional relevance, potentially affecting transmissibility and virulence. Experimental comparison of the effects of linoleic acid on SARS-CoV-2 variants, including emerging variants, is warranted.

ESCPE-1 mediates retrograde endosomal sorting of the SARS-CoV-2 host factor Neuropilin-1.

Endosomal sorting maintains cellular homeostasis by recycling transmembrane proteins and associated proteins and lipids (termed "cargoes") from the endosomal network to multiple subcellular destinations, including retrograde traffic to the trans-Golgi network (TGN). Viral and bacterial pathogens subvert retrograde trafficking machinery to facilitate infectivity. Here, we develop a proteomic screen to identify retrograde cargo proteins of the endosomal SNX-BAR sorting complex promoting exit 1 (ESCPE-1). Using this methodology, we identify Neuropilin-1 (NRP1), a recently characterized host factor for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, as a cargo directly bound and trafficked by ESCPE-1. ESCPE-1 mediates retrograde trafficking of engineered nanoparticles functionalized with the NRP1-interacting peptide of the SARS-CoV-2 spike (S) protein. CRISPR-Cas9 deletion of ESCPE-1 subunits reduces SARS-CoV-2 infection levels in cell culture. ESCPE-1 sorting of NRP1 may therefore play a role in the intracellular membrane trafficking of NRP1-interacting viruses such as SARS-CoV-2.

Structural insights in cell-type specific evolution of intra-host diversity by SARS-CoV-2.

As the global burden of SARS-CoV-2 infections escalates, so does the evolution of viral variants with increased transmissibility and pathology. In addition to this entrenched diversity, RNA viruses can also display genetic diversity within single infected hosts with co-existing viral variants evolving differently in distinct cell types. The BriSΔ variant, originally identified as a viral subpopulation from SARS-CoV-2 isolate hCoV-19/England/02/2020, comprises in the spike an eight amino-acid deletion encompassing a furin recognition motif and S1/S2 cleavage site. We elucidate the structure, function and molecular dynamics of this spike providing mechanistic insight into how the deletion correlates to viral cell tropism, ACE2 receptor binding and infectivity of this SARS-CoV-2 variant. Our results reveal long-range allosteric communication between functional domains that differ in the wild-type and the deletion variant and support a view of SARS-CoV-2 probing multiple evolutionary trajectories in distinct cell types within the same infected host.

Discovery of SARS-CoV-2 Mpro peptide inhibitors from modelling substrate and ligand binding.

The main protease (Mpro) of SARS-CoV-2 is central to viral maturation and is a promising drug target, but little is known about structural aspects of how it binds to its 11 natural cleavage sites. We used biophysical and crystallographic data and an array of biomolecular simulation techniques, including automated docking, molecular dynamics (MD) and interactive MD in virtual reality, QM/MM, and linear-scaling DFT, to investigate the molecular features underlying recognition of the natural Mpro substrates. We extensively analysed the subsite interactions of modelled 11-residue cleavage site peptides, crystallographic ligands, and docked COVID Moonshot-designed covalent inhibitors. Our modelling studies reveal remarkable consistency in the hydrogen bonding patterns of the natural Mpro substrates, particularly on the N-terminal side of the scissile bond. They highlight the critical role of interactions beyond the immediate active site in recognition and catalysis, in particular plasticity at the S2 site. Building on our initial Mpro-substrate models, we used predictive saturation variation scanning (PreSaVS) to design peptides with improved affinity. Non-denaturing mass spectrometry and other biophysical analyses confirm these new and effective 'peptibitors' inhibit Mpro competitively. Our combined results provide new insights and highlight opportunities for the development of Mpro inhibitors as anti-COVID-19 drugs.

The fatty acid site is coupled to functional motifs in the SARS-CoV-2 spike protein and modulates spike allosteric behaviour.

The SARS-CoV-2 spike protein is the first contact point between the SARS-CoV-2 virus and host cells and mediates membrane fusion. Recently, a fatty acid binding site was identified in the spike (Toelzer et al. Science 2020). The presence of linoleic acid at this site modulates binding of the spike to the human ACE2 receptor, stabilizing a locked conformation of the protein. Here, dynamical-nonequilibrium molecular dynamics simulations reveal that this fatty acid site is coupled to functionally relevant regions of the spike, some of them far from the fatty acid binding pocket. Removal of a ligand from the fatty acid binding site significantly affects the dynamics of distant, functionally important regions of the spike, including the receptor-binding motif, furin cleavage site and fusion-peptide-adjacent regions. Simulations of the D614G mutant show differences in behaviour between these clinical variants of the spike: the D614G mutant shows a significantly different conformational response for some structural motifs relevant for binding and fusion. The simulations identify structural networks through which changes at the fatty acid binding site are transmitted within the protein. These communication networks significantly involve positions that are prone to mutation, indicating that observed genetic variation in the spike may alter its response to linoleate binding and associated allosteric communication.

Molecular Simulations suggest Vitamins, Retinoids and Steroids as Ligands of the Free Fatty Acid Pocket of the SARS-CoV-2 Spike Protein*.

We investigate binding of linoleate and other potential ligands to the recently discovered fatty acid binding site in the SARS-CoV-2 spike protein, using docking and molecular dynamics simulations. Simulations suggest that linoleate and dexamethasone stabilize the locked spike conformation, thus reducing the opportunity for ACE2 interaction. In contrast, cholesterol may expose the receptor-binding domain by destabilizing the closed structure, preferentially binding to a different site in the hinge region of the open structure. We docked a library of FDA-approved drugs to the fatty acid site using an approach that reproduces the structure of the linoleate complex. Docking identifies steroids (including dexamethasone and vitamin D); retinoids (some known to be active in vitro, and vitamin A); and vitamin K as potential ligands that may stabilize the closed conformation. The SARS-CoV-2 spike fatty acid site may bind a diverse array of ligands, including dietary components, and therefore provides a promising target for therapeutics or prophylaxis.

Frontispiece: Molecular Simulations suggest Vitamins, Retinoids and Steroids as Ligands of the Free Fatty Acid Pocket of the SARS-CoV-2 Spike Protein

Drug Docking The binding of linoleate and other potential ligands to the recently discovered fatty acid binding site in the SARS-CoV-2 spike protein is investigated by James Spencer, Adrian J. Mulholland et?al. in their Research Article on page?7098.

A potential interaction between the SARS-CoV-2 spike protein and nicotinic acetylcholine receptors.

Changeux et al. (Changeux et al. C. R. Biol. 343:33-39.) recently suggested that the SARS-CoV-2 spike protein may interact with nicotinic acetylcholine receptors (nAChRs) and that such interactions may be involved in pathology and infectivity. This hypothesis is based on the fact that the SARS-CoV-2 spike protein contains a sequence motif similar to known nAChR antagonists. Here, we use molecular simulations of validated atomically detailed structures of nAChRs and of the spike to investigate the possible binding of the Y674-R685 region of the spike to nAChRs. We examine the binding of the Y674-R685 loop to three nAChRs, namely the human α4ß2 and α7 subtypes and the muscle-like αßγδ receptor from Tetronarce californica. Our results predict that Y674-R685 has affinity for nAChRs. The region of the spike responsible for binding contains a PRRA motif, a four-residue insertion not found in other SARS-like coronaviruses. The conformational behavior of the bound Y674-R685 is highly dependent on the receptor subtype; it adopts extended conformations in the α4ß2 and α7 complexes but is more compact when bound to the muscle-like receptor. In the α4ß2 and αßγδ complexes, the interaction of Y674-R685 with the receptors forces the loop C region to adopt an open conformation, similar to other known nAChR antagonists. In contrast, in the α7 complex, Y674-R685 penetrates deeply into the binding pocket in which it forms interactions with the residues lining the aromatic box, namely with TrpB, TyrC1, and TyrC2. Estimates of binding energy suggest that Y674-R685 forms stable complexes with all three nAChR subtypes. Analyses of simulations of the glycosylated spike show that the Y674-R685 region is accessible for binding. We suggest a potential binding orientation of the spike protein with nAChRs, in which they are in a nonparallel arrangement to one another.

Neuropilin-1 is a host factor for SARS-CoV-2 infection.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of coronavirus disease 2019 (COVID-19), uses the viral spike (S) protein for host cell attachment and entry. The host protease furin cleaves the full-length precursor S glycoprotein into two associated polypeptides: S1 and S2. Cleavage of S generates a polybasic Arg-Arg-Ala-Arg carboxyl-terminal sequence on S1, which conforms to a C-end rule (CendR) motif that binds to cell surface neuropilin-1 (NRP1) and NRP2 receptors. We used x-ray crystallography and biochemical approaches to show that the S1 CendR motif directly bound NRP1. Blocking this interaction by RNA interference or selective inhibitors reduced SARS-CoV-2 entry and infectivity in cell culture. NRP1 thus serves as a host factor for SARS-CoV-2 infection and may potentially provide a therapeutic target for COVID-19.

Free fatty acid binding pocket in the locked structure of SARS-CoV-2 spike protein.

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), represents a global crisis. Key to SARS-CoV-2 therapeutic development is unraveling the mechanisms that drive high infectivity, broad tissue tropism, and severe pathology. Our 2.85-angstrom cryo-electron microscopy structure of SARS-CoV-2 spike (S) glycoprotein reveals that the receptor binding domains tightly bind the essential free fatty acid linoleic acid (LA) in three composite binding pockets. A similar pocket also appears to be present in the highly pathogenic severe acute respiratory syndrome coronavirus (SARS-CoV) and Middle East respiratory syndrome coronavirus (MERS-CoV). LA binding stabilizes a locked S conformation, resulting in reduced angiotensin-converting enzyme 2 (ACE2) interaction in vitro. In human cells, LA supplementation synergizes with the COVID-19 drug remdesivir, suppressing SARS-CoV-2 replication. Our structure directly links LA and S, setting the stage for intervention strategies that target LA binding by SARS-CoV-2.